A proposal for a scalable universal bosonic simulator using individually trapped ions 
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We describe a possible architecture to implement a universal bosonic simulator (UBS) using 
trapped ions. Single ions are confined in individual traps, and their motional states represent the 
bosonic modes. Single-mode linear operators, nonlinear phase-shifts, and linear beam splitters can 
be realized by precisely controlling the trapping potentials. All the processes in a bosonic simulation, 
except the initialization and the readout, can be conducted beyond the Lamb-Dicke regime. Aspects 
of our proposal can also be applied to split adiabatically a pair of ions in a single trap. 



I. INTRODUCTION 



Coherently manipulated photons have been proposed 
to be a good candidate for testing the foundation of 
quantum mechanics [l|, performing quantum computa- 
tions 0, H[ , conducting high precision measurements [ij , 
and other many applications. However, because of poor 
sources, detection inefficiencies, and weak photon-photon 
interactions, implementing these proposals for large scale 
devices is very difficult. Some other well-controlled quan- 
tum system can simulate the photonic system if it also 
exhibits bosonic properties @. More specifically, a uni- 
versal bosonic simulator (UBS) should be able to repro- 
duce the evolution of a bosonic system under the most 
general form of Hamiltonian. This requirement is not 
too stringent as the evolution can be approximated to 
arbitrary accuracy by a sequence of basic operators that 
belong to a universal set @. Lloyd and Braunstein Q 
proved that the simplest universal set of basic operators 
comprised of all single mode linear operators, at least 
one multi-mode operator, and at least one nonlinear el- 
ement. Efficiently performing these basic operations is 
hence necessary for any implementation of the UBS. 

Ion traps are a suitable candidate for implementing a 
UBS, in which a high degree of controllability has been 
demonstrated. [8]. The motion of laser cooled ions is 
quantum in nature, and the excitations of the motional 
states, i.e. phonons, exhibit bosonic behaviour. The 
collective displacement and momentum of the ions are 
analogous to the quadratures of light fields. Any arbi- 
trary motional state can be created by combining tech- 
niques such as sideband transition M, parametric ampli- 
fication [Io| . and adiabatic passage [llj : in particular, the 
creation of Gaussian states [12[ and non-classical states 
[lH, [l3[ from the ground state have been experimentally 
demonstrated. When applied to non-ground states, some 
of these techniques can achieve single phonon linear or 
nonlinear operations. Interaction between phonon modes 
at the few-quanta level have also been observed. For 
example, nonlinear beam splitting on a single ion has 
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been performed by applying a Raman field [ljj ; coupling 
two phonon modes have also been demonstrated through 
the Coulomb interaction between two separately trapped 
ions [HI, EH ; or two ions in the same trap [Hi [HI ■ 

The idea of using trapped ions system for UBS was first 
proposed by Wineland et al. [l9(. Their simulator is a 
chain of ions trapped by a harmonic potential. If all three 
motional degrees of freedom are incorporated, a trap with 
N ions can at most offer 3iV phonon modes. This archi- 
tecture is in principle scalable, because the number of 
ions in a trap is not fundamentally limited. However, 
the addition of ions will narrow the frequency gap be- 
tween phonon modes; sideband transition should be con- 
ducted slowly to avoid significant errors [20( . In addition, 
a measurement on one phonon mode via resonance fluo- 
rescence may cause significant heating of the ion chain, 
which distorts the states of other phonon modes. These 
shortcomings limit the number of modes and the popu- 
lation of phonons that can be simulated accurately. 

The problem of an excess of ions in a single trap also 
appears in ion trap quantum computing [8j. Kielpin- 
ski, Monroe, and Wineland (KMW) [5l[ have proposed a 
modified approach in which ion qubits are stored in sep- 
arate locations of an array of traps, so that the manipu- 
lation on one qubit negligibly affects the others. Consid- 
erable advances in experimental realization of these ideas 
have been made in the past few years [22], [23J ; in partic- 
ular , entanglement gates have been performed on ions 
which were initially far-separated, and ions have been 
moved between traps. 

In this paper, we propose using the KMW architecture 
to implement a UBS on trapped ions system. We con- 
sider each ion to be stored in a separate harmonic trap, in 
which only one bosonic mode is present. All single mode 
operations can be conducted by either changing the stor- 
age trap potential or by laser manipulations. The linear 
beam splitter is based on the Coulomb interaction, which 
is the same principle as the kinetic energy exchange in 
Ref. [HI, [HI. However, we require the distance between 
the ions to be variable in order to speed up the process. 
Ions can be transported in specific trajectories that do 
not cause motional excitations (24|. The advantages of 
our scheme are that the quality of each operation is inde- 
pendent of the number of modes involved in the simula- 
tion, and the initialization and readout of any one mode 
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will not distort the others. 

We begin by presenting the setup and the physical 
model of our proposal in Sec. Ill Al Useful mathemati- 
cal techniques is introduced in Sec. Ill Bl The implemen- 
tation of single mode operations is then introduced in 
Sec. IIIII In Sec. IIVI we show that a linear beam splitter 
can be implemented by precisely combining and splitting 
two traps through changing the quadratic and quartic 
potentials. Initialization and readout of phonon states 
are presented in Sec. |V] This article is concluded in Sec. 
[VT]with some discussions. 



II. IONS IN HARMONIC TRAPS 
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FIG. 1: The configuration of our ion trap UBS is an array of 
storage traps. Only one ion is trapped in each trap, which 
is a harmonic well. The traps are separated by a distance L, 
which is large enough to prevent disruption from the others. 
The position displacement of the ith ion, Xi, is accounted with 
respect to the trap centre. 



A. Model 

We assume ions are tightly trapped in the y and the z 
directions by a strong ac held while a weaker dc potential 
is applied along the x direction. We assume this config- 
uration would effectively restrict the ion to move along 
only the x direction because the excitations in other di- 
rections are negligible. 

The configuration of our system is schematically shown 
in Fig. Q] Ions are trapped in an array of harmonic 
storage traps, and only one ion is present in each trap. 
The distance between the equilibrium position of two 
neighbouring traps is L, which is sufficiently large that 
Coulomb coupling between the ions can be neglected. 
Hence the total Hamiltonian of the system is given by 
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where the subscript n denotes the quantities belonging 
to the ion in the nth trap; x n is the operator of the nth 
ion's location measured from the equilibrium position of 
the nth trap. The annihilation and the creation operators 
of the phonon mode of the nth ion are given by 



2h 



x n + i 



1 



2hmujQ 



/ muiQ „ 
./ 2h x„ 



1 



2hmujQ 



Pn 



Pn 



(2) 



The ions are cooled to both the motional and electronic 
ground state before an input state is created. The trap- 
ping potentials will be varied, but the potentials should 
return to the original form in Eq. ([TJ after each opera- 
tion. Thus, an operation is characterized by the transfor- 
mation of the motional state in the interaction picture, 
i.e. 



\Mn = s\mo)) 



(3) 



where S = exp(iHot/h)Us is the S-matrix; Us is the 
evolution operator in the Schrodinger picture. Then the 



annihilation operator of a phonon mode in the interaction 
picture is transformed as 



In -> S^a n S = Ula n Use 



iu) T 



(4) 



We will omit the n in future discussions of single mode 
operation. 



B. Lewis- Riesenfeld theory 

The Schrodinger equation with a time dependent har- 
monic potential will frequently be encountered in future 
sections, i.e. 



id t \^) = H(tm 
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Lewis and Riesenfeld analyzed the problem of time vary- 
ing harmonic oscillators by considering the invariant op- 
erator given by the following expression [2o1 . [26| , 



±{t) 



(bp — mbx) 2 1 



2m 



1 



+ -mu> 2 x 2 = hu>o ^A t (f)A(t) + 2 j ■ 

where a dot denotes a time derivative. The dimensionless 
function b(t) satisfies the auxiliary equation 



b + u; 2 (t)b-^-=0 
b 6 



(7) 



and the operators A(t) and A^(i) are the raising and low- 
ering operators of the eigenstates, |A„,i), of i.e. 



±{t)\X n ,t) = A„|A„,t) , 
A(t)\X n ,t) = Vn\X n -!,t) 
At(i)|A n _i,i) = V^\K,t) , 



(8) 
(9) 
(10) 



with A n being the corresponding eigenvalues. The invari- 
ant operator is so called because [I(t),H(t)] = 0. Hence 
the values of X n remain unchanged, and \X n ,t) are the 
solutions of Eq. ([5]) during the evolution, i.e. 



(X m ,t\H(t)\X n ,t) = , if n^m 



(11) 
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The evolution operator of the system from time t to t' is 
given by the expression 

oo 

&(t,o = E ei(n+ * )(e(t) " e(t ' )) i A «'*)< A «'* / i ' ( 12 ) 



where the phase is given by 



Q(t) 
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When the harmonic well is static, i.e. lu is a constant, 
the general real solution of Eq. ([7]) is 



- sj cosh 6 + sinh S sin(2cj£ + ip) , (14) 

where d and tp are constant parameters. In our case, we 
are only interested in the operations where the trapping 
potential are steady before and after the operation, i.e. 



u(ti) = u(tf) = oj , 



(15) 



where ti and 1 1 are the starting and ending time of the 
operation. For simplicity, we assume both 6 = and ip = 
at the beginning, such that b{tj) = 1. In general, the 
values of 5 and tp have to be determined by integrating 
Eq. ([7]). The invariant operator X(t) is identical to H (t) 
at t = ti, so we have 

k(ti) = a . (16) 
After the operation, the lowering operator becomes p5| 



where 



k(t f )=r,(t f )a + t(t f )tf 
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The absolute magnitude of r)(t) and are 

5 S 
\r,(t)\ =cosh- ; |C(t)| =sinh- , 
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(17) 

(18) 
(19) 

(20) 
1 is sat- 



which the normalization condition |?7| 2 
isfied automatically. 

The action of the harmonic potential variation is ac- 
counted by the evolution operator, i.e. O = U(tf,ti). It 
is readily seen that the raising operator is related to the 
annihilation operator as 

k(t f ) = e^ e ^- e ^U(t f ,U)aUHt f ,U) . (21) 

Using Eq. (j4]) , we deduce the transformation of annihila- 
tion operator as 

-^(t / ) e - i ( e ( t /)- e ( t «)-"o(t / -t i )) a t . (22) 

A generalization of Lewis-Riesenfeld theory, including 
motion of the trap centre, is given in Appendix IA1 



(a) 



(b) 
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FIG. 2: (a) A displacement operator is implemented by chang- 
ing the trapping centre of the harmonic well, (b) A phase-shift 
operator or a squeezing operator is implemented by varying 
the harmonic potential strength, (c) An extra quartic poten- 
tial is applied to implement the nonlinear phase gate. 



III. SINGLE MODE OPERATIONS 

Any single mode linear operator can be achieved by al- 
ternatively applying the displacement operators, phase- 
shift operators, and squeezing operators @, H3- For ion 
trap bosonic simulations, each of these operators could be 
implemented by constructing specific Hamiltonians using 
laser interaction. However, the accuracy and speed are 
limited by the validity of the Lamb-Dicke approximation 
(LDA), unless complicated higher order corrections are 
considered 0, EH ■ We consider an alternative ap- 
proach that the operators are implemented by varying 
the harmonic trapping potentials. In addition, applying 
perturbatively a quartic potential to the storage trap can 
achieve a nonlinear phase gate, which is a nonlinear op- 
erator that comprise the universal set of operators. Both 
harmonic potential and quartic potential can be implc- 
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mcnted in experiments 
assumed to operate from t 
the operations is shown in Fig. [21 



All of the operations are 
Q to t = T. Summary of 



A. Displacement Operator 

A displacement operator, T>{a) = exp(aa^ — a* a), 
transforms the annihilation operator as 



a -> (a)aV(a) =a + 



a 



(23) 



where y // 2h/mLUoRe(a) is the shift of the ion's position; 
y/2hmujolm(a) is the shift of the ion's momentum. The 
operator can be achieved by applying two radiation fields 
with close frequencies, to induce a Raman transition. If 
the frequency difference of the Raman fields is resonant 
to the red sideband of the ground electronic state, the 
effective Hamiltonian is then proportional to a + at 
the first order of the LDA and the rotating wave approx- 
imation (RWA) , while the internal state remains at the 
ground state [lj, [H, [2(| ■ 

Another way of performing the displacement operator 
is to move the harmonic trap, i.e. replacing the storage 
trap Hamiltonian by the displacement operator Hamilto- 
nian 



p 2 1 



(24) 



where s(t) specifies the path of the trap centre . Assume 
the trap centre initially located at the origin, i.e. s(0) = 
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0, we require the trapping centre returns to the origin 
after the operation, i.e. s(T) — 0. After the operation, 
every coherent state, \x), will become 



Ix) 



muJo 
~2h 



sit) exp(iu)ot)dt 



-iu T 



(25) 

up to a global phase that will not affect the simulation 
result j29|. Since the result applies to every coherent 
state, the annihilation operator transforms according to 
Eq. (0} is 



mujQ 
2h 



s(t) exp(iujQt)dt = a + a[s(t)] , (26) 



where the displacement a is a functional of s(t). We 
note that many different s(t) can produce the same a. A 
method to obtain an s(t) for a specific a is described in 
Appendix [B] 



B. Phase-Shift Operator 

A phase-shift operator P((j>) = exp(—i<pa^d) trans- 
forms the annihilation operators as 



a -S- ^((fyaP^) = de 



(27) 



The operator can be implemented by applying a laser 
field to the ion, where the laser frequency is detuned from 
any electronic transition. The ac Stark effect produced 
by the field will disturb the electronic ground state en- 
ergy, and hence induces a phase shift. 

Wineland et at. [3] implemented the phase-shift oper- 
ator by applying a perturbative quadratic potential, ex . 
Under the RWA, the term heiaa^ 4- a) a) /2mu>o domi- 
nates the perturbation and produces a phase shift. In 
this approach, the validity of the RWA requires that the 
strength of the perturbation, he/2mujQ, be much smaller 
than the energy of a phonon hu>o- The requirement im- 
plies that the duration of a phase-shift operator, which 
scale as \Jm/e, should be much longer than 1/wo- 

The operation time can be reduced to the same or- 
der as l/wo if the strength of the harmonic trap changes 
non-perturbatively, i.e. replacing the storage trap Hamil- 
tonian by the phase-shift operator Hamiltonian 



~2 



f 



H P = ^- + ^-muj 2 (t)x 2 
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(28) 



The trapping frequency should be the same as that of 
the storage trap before and after the operation, i.e. 
w(0) = u{T) = ojq. The effect of this potential can be 
considered analytically using the Lewis-Riesenfeld theory 
fi"J introduced above. Because the amplitude of each 
motional Fock state remains the same after the phase- 
shift operator, it implies that |A ra ,0) = \X n ,t >T) = \n) 
up to a phase, where \X n ,t) is an eigenstate of the in- 
variant operator at time t. This criterion is equivalent 



to require the corresponding auxiliary function of u>(t), 
&^(t), which follows Eq. ([7]), to satisfy the boundary con- 
ditions 



M* < o) = 1 ; hit > r) = 1 



(29) 



Then the final lowering operator will consist of only the 
annihilation operator, i.e. r)(t > T) = 1, ((t > T) = 0. 
According to Eq. (|2"2"|) , the overall effect of the harmonic 
potential variation will transform the annihilation oper- 
ator as 



ae 



He(T)+uj T) 



(30) 



which is obviously a phase-shift operator. 

There is no unique form of ui(t) and b^t) that satisfies 
all the above conditions, therefore we have the freedom 
to choose u)(t) in a manner that is convenient in practice. 
Alternatively, we can initially guess a b(t) and obtain the 
corresponding uj(t) for the experiment. A possible choice 
is b<t,(t) = 1 - fcexp(-(i - T/2) 2 /a 2 ) , where l/a < T 
is the characteristics time scale of the operation; k is 
chosen to produce the desired phase shift. There is no 
fundamental limitation on the magnitude of cr, so our 
phase-shift operator implementation can be processed in- 
definitely fast, which can even be faster than 1/wo if the 
apparatus permits. 



C. Squeezing Operator 

A squeezing operator S(g) = exp ((g*a 2 — ga j ' 2 )/2) 
transforms the annihilation operator as 



a —> (g)aS(g) = cosh \g\d — sinh \g\a? 

\g\ 



(31) 



A squeezing operator is usually implemented by applying 
an interaction that the Hamiltonian involves second or- 
der terms in the annihilation and the creation operators, 
i.e. a 2 and cl* 2 . Such a Hamiltonian can be realized by 
Raman interaction with effective frequency 2loq, which is 
initiated by applying two radiation fields with 2cu differ- 
ence in frequency. However, the magnitude of the poten- 
tial generated by the Raman fields must be much smaller 
than 2haJo to satisfy the RWA. The operation time of the 
squeezing operator will then be much longer than 1/wo- 
Our approach is to use a time varying trapping po- 
tential, i.e. replacing the storage Hamiltonian by the 
squeezing operator Hamiltonian 



P 



1 



2m 2 



(32) 



The trapping frequency is required to return to that of 
the storage trap after the operation, i.e. w(0) = uj{T) = 
ujq. The operation will transform the annihilation opera- 
tor according to Eq. (|22[) , which is apparently a squeezing 
operator if |£(T)| 7^ 0. The magnitude of the squeezing 
parameter is given by \g\ = 8/2 according to Eq. (|20l) . 
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which the S has to be obtained by integrating the auxil- 
iary equation Eq. with the ui(t) applied in the oper- 
ation. The phase of the squeezing operator, g/\g\, comes 
from the complex nature of £(T) and r)(T), and the phase 
shift of (8(T) + uj a T) in Eq. @. The effect of this op- 
eration can be accounted by analytically solving the time 
dependent harmonic oscillator, so the operation time is 
not limited by the validity condition of the RWA. 

The applied potentials of H$ and Hp are both time 
varying harmonic wells; the only difference is the time 
variation of the trapping frequency, u){t). Unless u>{t) is 
specially designed as shown in Sec. IIII B[ the operation of 
varying the harmonic potential would be a squeezing op- 
erator with some squeezing parameter crucially depend- 
ing on uj(t). The uj(t) that generates a specific squeezing 
parameter can be obtained by the method detailed in 
Appendix [Cl 



D. Nonlinear Operator 

Nonlinear operators transform an annihilation opera- 
tor to an operator involving quadratic and higher order 
terms in a and . It can be achieved by applying a 
Hamiltonian that is at least third order of a and . One 
implementation is to exert a radiation field that is reso- 
nant to high order sideband frequencies. For example, a 
Raman field with the effective frequency 3wo can provide 
a Hamiltonian scaling as (a 3 — a^ 3 ) at the third order ex- 
pansion of the Lamb-Dicke parameters [12j . The major 
problem of this approach is that the validity of both the 
LDA and the RWA have to be satisfied, so the undesired 
Hamiltonian are suppressed. The consequence is that the 
power of the radiation field is constrained, which limits 
the speed of the operation. Nonetheless, our architec- 
ture can facilitate this laser-mediated nonlinear operator, 
because the mode spectrum is simplified since only one 
bosonic mode is exhibited in each trap. 

A nonlinear operator can also be implemented by 
switching on an additional quartic potential, i.e. 



H 4 (t) = H + Vi(t) = H + T(t)x 4 



(33) 



In the interaction picture with respect to H$, the quartic 
potential becomes 



Vl 



T{t)h 2 



(6a t2 a 2 + 12a) a + 3) 



(2towo) 2 

+a) 2 {Aa)a + 6)e i2 <"°' + a t4 e l4 " ot + h.c.l .(34) 



If the variation of F(t) is slow enough, the off-resonant 
terms can be eliminated by the RWA; the only effective 
terms are 

H N ee L {t)H " {&a) 2 a 2 + 12&t& + 3) . (35) 

(ZTOWo) 

Applying the quartic potential from t = to T, the 
S-matrix of the operation in the Schrodinger picture is 



given by 



where 



S 4 



— e -i/j(T)(6a t2 a 2 + 12a+a) 



JO 



(2mwo) 2 



di' 



(36) 



(37) 



We have neglected the unimportant global phase, and 
have employed the fact that [H ,Hn]. 

The speed of the operation is mainly determined by the 
validity of the RWA. According to Ref. [3(|, applying the 
RWA is to collect the leading order terms in a series ex- 
pansion of time-averaged Hamiltonians. The sufficient 
condition for a valid series expansion is that the largest 
eigenvalue of Hm / H should be much smaller than the off- 
resonant frequencies, which are multiples of loq in our 
case. Although / h has unbounded eigenvalues, the se- 
ries expansion is still valid if the maximum phonon num- 
ber, rt m axj m each mode is small. To estimate the RWA 
validity condition, we approximate J 7 (t)fi 2 /(2mwo) 2 by 
Hfi(T)/T because V4 is slowly varying. The maximum 
eigenvalue of Hn in our simulation is hence n% aax IJ<{T)/T, 
which gives a valid RWA when 



UJqT 



< 1 



(38) 



Arbitrary evolution of a single phonon mode can be 
realized by repeatedly applying the linear operators and 
the nonlinear phase-shift operator Q- The idea is to 
consider the intertwinement of the infinitesimal evolution 
of two Hamiltonian y and 2, i.e. 

e iyst/H e iZ6t/H e -iyst/H e -iZ6t/h _ e i[y,z]St 2 /h 2 +0 ^3\ 

(39) 

If the 0(St 3 ) terms are neglected, the resultant opera- 
tion of this sequence is effectively the evolution of a new 
Hamiltonian [3^,2]. Evolution of new classes of Hamil- 
tonian will be generated by using this trick again; the 
evolution of the desired Hamiltonian is eventually pro- 
duced. 

The corresponding Hamiltonian of a linear operator 
consists of second or lower orders of a and a'. For our 
selection of linear operators, the Hamiltonian of T*(a), 
V{(j>)i S(g) are respectively proportional to a, (r a, a 2 and 
their Hermitian conjugates. Hamiltonians with higher 
than second order of a and cannot be generated by 
intertwining the linear operators, therefore the nonlinear 
phase-shift operator is required to implement UBS. 



IV. TWO-MODE OPERATION 

In Sec. IIII1 we examined the techniques to realize single 
mode operations. We now turn our attention to mixing 
two modes, analogous to a beam-splitter in optics. A 
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FIG. 3: Displacement of ions and variations of potentials dur- 
ing a beam splitter operation. The origin is denned as the 
mid-point between two storage traps. Step I, ions are trans- 
ported by harmonic well from storage traps to pick-up posi- 
tion. Step II, double well is switched on to pick-up ions. Step 
III, separation of double well shrinks. Step IV, the separation 
of double well extends and brings the ions back to pick-up po- 
sition. Step V, double well is switched off, and the harmonic 
wells pick-up ions and bring them to storage traps. 



position and momentum operator; Xi and Pi = mXi 
are the classical position and momentum of the ith ion; 
qi and iti are the operators accounting for the quantum 
fluctuation of position and momentum about Xi and Pi 
respectively. The aim of our beam splitter is to transform 
the quantum fluctuations of the two ions according to 
Eq. (|40|. whereas the ions will be classically stationary 
at the storage traps before and after the operation, i.e. 

Xi(~) = = ~\ 5 X 2 (~) = X 2 (|) = | ; (42) 



Px(- 



2' 



- T 



P 2 (- 



1< 

2' 



- T 



. (43) 



At the storage traps, the quantum fluctuations are the 
same as the phonon quadratures defined in Eq. ([T]), i.e. 



storage 



I storage 



Pi 



(44) 



The core component of our phonon beam splitter is 
a double well potential with varying well separation. It 
can be constructed by a quartic potential, A(t)X 4 , and a 
harmonic potential, B(t)X 2 , which can be implemented 
in experiments [28]. The evolution of the motional state 
follows the Schrodinger equation ihdt\ip) — HBsit)^), 
where 



linear beam splitting operation, B(9, 
boson modes as 



transforms two 
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(40) 



A 50:50 beam splitter corresponds to the case 9 = ir/2. 

The whole beam splitter process is summarized 
schematically in Fig. [3] Displacing harmonic wells are 
applied from t = —T/2 to transport the ions from the 
storage traps to the pick-up distance and then switched 
off at t = —T'/2, a double well potential is immediately 
switched on to relay the transportation. The separation 
of the double well shrinks and then expands; the two ions 
are brought to proximity and then separated. The two 
encoded phonon modes interact via the Coulomb inter- 
action between the ions. The double well potential fi- 
nally separates the ions to the pick-up distance and then 
switched off at t = T'/2, two moving harmonic wells are 
switched on immediately to transport the ions back to 
the storage traps at t = T/2. 

In this section we set the origin of position, X = 0, 
is the mid-point between the two storage traps. We as- 
sume the system is both spatially and dynamically sym- 
metric about the origin. To simplify the discussion, we 
will separate the evolution of the classical motion and the 
quantum fluctuations by defining 



Xi — Xi 



Pi = Pi 



(41) 



where the subscripts £ = 1,2 denote the ions involved 
in the beam splitter operation; X and P are the total 



P 



p: 



H BS (t) = -± + -± + B(t)(Xt + X. 



2m 2m 
-A(t)(Xt 



ATre Q (X 2 -X 1 ) 



;(45) 



is the total wave function. We have assumed the ions 
cannot tunnel pass each other due to the strong Coulomb 
repulsion, so X 2 > X\ is always true. 

In terms of the variables in Eq. (|4ip . the Schrodinger 
equation becomes 



(46) 



where "Ho and "Hi collect the terms with the zero and 
first order quantum fluctuations, and Hb contains the 
rest; \tp) is the state of the quantum fluctuations. The 
first term H = Pf /2m +_P 2 2 /2m + A(t)(Xf + Xf) + 
B(t)(Xl + Xl) + e 2 /47re (A 2 - X x ) is the total mechani- 
cal energy of the system; it only contributes to an unim- 
portant global phase. The second term Ti.\ vanishes if 
the classical equations of motion are satisfied, i.e. 



Xi = Pi/m , 

A = -4A(t)Xf - 2B(t)Xi 



e 2 X : 



4 7 re |A J |(A 2 -A 1 ) 2 



Because of the symmetry, we have Xi(t) = —X 2 (t) and 
Pi(t) = —P 2 [t). The classical separation between the 
ions is defined as r — X 2 — X\, then the equation of 
motion reduces to the following: 



m 



2B(t) 



Aneomr 2 



(47) 
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If the quantum fluctuation of position is much smaller 
than the separation of ions, i.e. yj (5q 2 )/r <C 1, then Hb 
can be approximated by a quadratic Hamiltonian, viz 



e 2 (Q2 - <7i) 2 



(48) 



Instead of staying in the bases of individual ions, it 
is advantageous to consider the centre-of-mass mode (+ 
mode) and breathing mode (- mode). The position and 
momentum operators are defined as 



q 2 ±qi 7T2 ± 7Tl 

q± = — — ; n ± = 



V2 



V2 



(49) 



In the new basis, Eq. (|48[) can be decoupled as two har- 
monic oscillators, 

% = ^+^(^+g+^(t)e, (50) 

where the mode frequencies are 



3A(t) , 2B(t) 



;_(t) = + 



(51) 
(52) 



27reomr 3 

The annihilation operators of the modes are defined as 



. muj± „ 



i 



2hmuj± 



(53) 



Changing the magnitude of the double well will result 
in either squeezing or phase-shifting the + mode and the 
- mode. Because there is no excitation after a beam split- 
ting, the double well operation should give only a phase 
shift. For simplicity, we assume the quartic and the har- 
monic potentials are adjusted to produce a constant w+, 
i.e. 



3A(t) 2 



r\t) 



2B(t) 



(54) 



According to Eq. (|2"2")l . the + mode remains unchanged 
after the operation, i.e. A+ — > A+. 

We require the double well operation is merely a phase- 
shift operator on the breathing mode, i.e. the anni- 
hilation operator transforms as A- —> e~ l6 A-, where 
9 = (-e(T) - ujqT) according to Eq. (J25J. The phonon 
modes of individual ions then transform as 

1 6 

Ai —j={A+ - e~ l6 A-) = cos -A + isin-„4 2 ; (55) 

V 2 2 2 

A 2 -> -^(-4+ + e~'' e .4_) = i sin - A + cos - A , (56) 

V 2 2 2 



which is a beam splitter operation. We have neglected 
the unimportant global phase e~ l9 ^ 2 , and the phase <f> 
in Eq. (j40|) can be rectified by applying local phase-shift 
operators. 

The pick-up distance should be large enough that u;_ 
and w+ should be roughly the same according to Eq. (|52|) , 
i.e. LJ-(T'/2) =w_(T'/2) = ojq. To produce only a phase 
shift on - mode, the should produce an auxiliary 

function b(t) that satisfies b(-T'/2) = b(T'/2) = 1. We 
can control the strength of the quartic and harmonic po- 
tentials, A(t) and B(t), to produce an appropriate W-(t) 
while preserving u>+(t). The appropriate time variations 
of A(t) and B(t) exist and are not unique; they can be 
chosen in a manner that is convenient in practice. In Ap- 
pendix |D1 we suggest a method to obtain A(t) and B(t) 
from a speculative b(t). 

Before and after the double well operation, the ions are 
transported back and forth between the storage traps and 
the pick-up distance (step I and V). If both the transport- 
ing harmonic potentials and the double well potential can 
be switched on and off quickly, the pick-up process can 
be conducted smoothly that the phonon states will not 
be disturbed. The pick-up distance is arbitrary; it can 
be chosen in a manner that is favourable in experiments. 

The ions' classical velocity in the double well operation 
is determined by the choice of A(t) and B(t). The veloc- 
ity at the pick-up distance, X,(-T'/2) and X l (T'/2), is 
obtained by integrating the equation of motion Eq. (|47[) 
and the condition of a vanishing velocity at the turn- 
ing point, i.e. r(0) = 0. In step I, the transporting 
harmonic wells should increase the classical velocity of 
the ions from at the storage trap to Xi(—T'/2) at the 
pick-up distance. Similarly in step V, the transporting 
harmonic wells should reduce the classical velocity from 
Xi(T'/2) at the pick-up distance to at the storage traps. 
The classical position and momentum of the ions in the 
transportation process depend on the trap centre, 
of the transporting harmonic well, which the Hamilto- 
nian is H T = Pf/2m + mu} 2 (Xi - &(t)) 2 /2. Using the 
variables in Eq. (14"TT) . the classical equation of motion is 



mXi(t) 



"1 (&(*) 



(57) 



where the exact solution can be found in Ref. [29j. Ap- 
propriate £i(t), which produces the Xi and Pi that match 
the boundary conditions at t = -T/2,-T'/2,T'/2,T/2, 
can be obtained by the inverse-engineering method pre- 
sented in Ref. [2J|. The evolution of the quantum fluc- 
tuations in the transportation follows the Hamiltonian 



ij ^ ^ 1 - 
Ut= 2^ + 2 mUQq 



2„2 



(58) 



Obviously, the phonon states will not be disturbed by Ht 
and hence the transportation. All in all, the operation 
from step I to step V realizes a phonon beam splitter, i.e. 
Eq. (|40[) . on the phonon modes in neighbouring storage 
traps. 
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Now let us consider the errors in the beam splitter op- 
eration. The transportation in step I and V would cause 
small error if the harmonic well is sufficiently accurate. 
The pick-up process is assumed to be fast enough that 
does not cause significant error. Thus, most of the er- 
ror is expected to come from the double well process. 
The main problem is that the anharmonic terms in Kb, 
which are the terms with at least third order of quan- 
tum fluctuations and are not covered by H2, would pro- 
duce unwanted excitation or phase shift. The error is 
expected to be magnified if more phonons are involved in 
the beam splitter, because the significance of the quan- 
tum fluctuation is characterized by \J (q 2 )/r, which in- 
creases as phonon number. A faster operation also causes 
more significant errors, because some of the anharmonic 
terms would be suppressed by RWA if the operation is 
slow. The evolution of phonon states is calculated in Ap- 
pendix [E] when some anharmonic terms are included in 
the Hamiltonian. We find that a beam splitter can be 
conducted as fast as 3 /is if there are less than 8 phonons 
involved in each Ca + ion trap. 

The technique of double well operation can also be 
applied to split adiabatically a pair of ions in a single 
trap, details are provided in Appendix IF1 



V. INITIALIZATION AND READOUT 

Manipulation of trapping potentials is sufficient for ini- 
tializing arbitrary Gaussian phonon states, which can 
be created by applying linear operators to the ground 
motional states. However, non-Gaussian phonon states, 
such as Fock states and the Schrodinger's cat state, have 
to be generated by laser interaction. When comparing 
with the previous UBS proposal that incorporates a chain 
of ions in a single trap [19L l20j , our architecture simplifies 
the laser-mediated state initialization. Because the ions 
are separately trapped, sophisticated techniques to pre- 
vent the laser operation from disturbing other ions, such 
as using composite pulses and shielding, are no longer 
necessary. 

Information about the phonon states can be extracted 
by the three measurement schemes suggested in Ref. 
20]: using adiabatic transfer, post-selection techniques, 
and multiple electronic states. The adiabatic transfer is 
achieved by exerting a Raman pulse whose frequency is 
slowly increased from lower to higher than a sideband fre- 
quency [III ■ The phonon-dressed electronic states trans- 
fer as \g n) |e (n. — 1)) while \g 0) remains unchanged. 
The fluorescence measurement scheme in ion trap quan- 
tum computing is conducted after one round of adiabatic 
transfer, i.e. by applying a strong laser pulse resonant to 
the transition frequency between \g) and some other un- 
stable electronic state of the ion; significant fluorescence 
is detected if the electron is in the state \g). A positive 
measurement outcome is produced if the original phonon 
state is |0). This procedure is equivalent to the non- 
distinguishing number detector in optical experiments. 



If the frequency detuning and the Rabi frequency of the 
Raman field are both tuned, an adiabatic transfer can be 
achieved for a single mode with 0.99 fidelity in as fast as 

80 Ais [Hill]. 

If more rounds of adiabatic transfer and car- 
rier pulses are applied, a fluorescence measurement 
can produce the projection- value measure (PVM), 
U2n=o \ n )( n \,J27= m +i l n )Hi' wni ch can distinguish if 
the phonon population is more or less than an integer 
m. If there is another meta-stable electronic state, then 
any Fock state |m) can be deterministically measured by 
the PVM {|m)(m|,I — |m)(m|}, where I is the identity 
matrix. Incorporating with the displacement operator, 
any coherent state |a) can also be deterministically dis- 
tinguished by the PVM {|a)(a|,I - |a)(a|} [Hj]. To the 
best of our knowledge, detectors with these PVM have 
not yet been developed in optical experiment. 

The post-selection measurement method is applicable 
only if a negative fluorescence measurement causes neg- 
ligible distortion on the resultant phonon state. Post- 
sclcction quantum information protocols, such as the lin- 
ear optics entanglement gate [2| , can be achieved by using 
this measurement method. A sideband pulse is applied 
to transit some unwanted superposition to |e), while the 
desired outcome state remains in \g). The unwanted su- 
perpositions will be removed if the fluorescence measure- 
ment with respect to |e) gives a negative result. Repeat- 
ing this process will increase the amplitude of the desired 
outcome state in the residual phonon state. 

If the maximum number of phonons in a mode is known 
as ri max , and there are n max meta-stable electronic states 
available for manipulation, each Fock state with phonon 
number smaller than or equal to n max can be projectively 
measured, i.e. phonon number distinguishing measure- 
ment. The principle is to associate each phonon Fock 
state to an electronic state; the fluorescence measure- 
ment is then conducted on the electronic states one by 
one. The sequence of pulses employed to measure a mo- 
tional state with n max = 2 is given in [201 ] . 

In all the three measurement schemes, our architecture 
is more favourable than the previous proposals in Ref. 
(l9l [20] ] that traps a chain of ions in a single harmonic 
potential. Because the ions are individually trapped, the 
recoil of an ion after fluorescence measurements do not 
distort other phonon modes. The spectral distribution 
of resonance is also simplified because only one mode is 
present in each trap; the speed of sideband transition is 
hence increased as a stronger pulse can be used without 
the concern of accidental mode mixing. 



VI. CONCLUSION 

We have described a possible architecture to imple- 
ment the universal bosonic simulator using individually 
trapped ions. The excitation of an ion's quantized mo- 
tion can simulate a bosonic mode. Linear single mode 
operations can be realized by changing the strength and 
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the trapping centre of the harmonic potential at the stor- 
age trap. Nonlinear phase-shift operator can be imple- 
mented by perturbatively exerting a quartic potential. 
Linear beam splitter is implemented by a double well 
potential with varying separation; the Coulomb interac- 
tion between the ions ensues from the interaction between 
phonons. By applying the operators alternatively, arbi- 
trary bosonic evolution can be effectively simulated @, @] ■ 
All linear and nonlinear operators can be conducted 
without laser interaction, hence the speed of the opera- 
tions are not limited by the Lamb-Dicke approximation. 
However, all the laser-mediated techniques of the previ- 
ous proposals in Ref . [lij , which employ a chain of ion 
in a single trap, are applicable to our architecture. Be- 
cause only one phonon mode is associated to each trap, 
the resonant frequency spectrum is much simpler than a 
trap with multiple ions. The requirement of the rotating 
wave approximation is also less stringent due to the ab- 
sence of accidental mode mixing. In addition, measuring 
a phonon mode by laser keeps the other phonon modes 
undistorted. 

Provided that the quality and the controllability of the 
harmonic potential are fine enough, the operation time 
of the single mode linear operators has no fundamen- 
tal limit, and a two- mode operator can be implemented 
within tens of l/wo if the phonon number in each inter- 
acting mode is less than 10. Although the speed of the 
nonlinear operator is relatively slower, various interesting 
bosonic phenomena can be investigated using only linear 
operators. 

Recently, Aaronson and Arkhipov [33[ proposed that 
if there exists a classical algorithm that efficiently sam- 
ples the probability distribution of a linear bosonic net- 
work, then the polynomial hierarchy would collapse to 
the third level, which is generally believed to be impos- 
sible in computer science [34j]. In other words, an ap- 
proximate bosonic sampler should be a machine having 
post- classical computing power; building such a machine 
can verify the power of quantum computer. Aaronson 
and Arkhipov suggest that a boson sampler involving 
n = 10 to 50 bosons and about ~ n 2 to ~ n 5 logn modes 
suffice to achieve the goal. The scalability, speed, the- 
oretical quality, and measurement flexibility of our ion 
trap bosonic simulator architecture offers the possibility 
of demonstrating such post-classical computation in the 
not too distant future. 

We thank Marcos Villagra, Christian Weedbrook, and 
Andrew White for useful comments. The authors would 
like to acknowledge support from the NSERC CREATE 
Training Program in Nanoscience and Nanotechnology. 



Appendix A: Generalized solution of harmonic 
oscillator 

We have introduced the method to individually imple- 
ment the displacement operator and the squeezing oper- 
ator. In fact, both operators can be implemented in one 



operation, which employs a harmonic trap with both the 
strength and the centre of the trap are varying, i.e. 

md t \m) = + \mu 2 {t){x - s{t)A \m) ; (ai) 

where \ip(t)) is the solution of the equation. This gen- 
eralized harmonic oscillator has been investigated in the 
context of the evolution of number states [35| . and the 
variation of quadratic operators in the Heisenberg pic- 
ture [36l |. We here formulate the previous results to fit 
in our purpose of bosonic simulation, i.e. to consider the 
evolution of the annihilation operator in the interaction 
picture according to Eq. (|4]). 
We define the displaced state as 

\ X {t))=&{(5{t)e-™ ot )\m) > (A2) 

where T> is the displacement operator; the displacement 
is given by the expression 

me -^ = ^ M+i ^Z Mt) . (A3 ) 

If x c and p c are respectively the classical position and 
momentum of the ion, i.e. they obey the classical equa- 
tion of motion 

± c (t)=p c (t)/m ; p c (t) = -mcu 2 (t)(x c (t) - s(t)) , (A4) 

then the displaced state follows the equation 

m\x(t)) = + \^{t)sA \x(t)) , (A5) 

which can be solved by the Lewis-Riesenfeld theory. An 
unimportant global phase has been neglected in the above 
derivation. 

The operation runs from t = to T. The trap centre 
and strength should return to that of the storage trap at 
t = T. The evolution operator of the displaced state in 
the Schrodinger picture, U x .s, transforms the annihila- 
tion operator as 

Uls^s = V*(T)e^a - C(T)e-^a^ . (A6) 
The evolution operator of \ip(t)) is related to U x ,s by 

Us(t) = V(J3(t)e- iu ° t )U x ,s(t). (A7) 

Using Eq. (Q} , we can deduce the transformation of the 
annihilation operator in the interaction picture is 

a -> ry*(T)e 4 ( e ( T '+^ T )a 

_£( T ) e -i(e(T)-a, T) a t + /3( T ) . (A8) 

If b(t) is known, /3(T) can be obtained as (3|| 

P{T) = l \f^o Jo mU32{tHt)e ^ e{t)dt 

x (V(T)e l ( e ( T >+"« T > + ((T)e- l( - e ^'^ T ^ (A9) 
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The total operation is obviously a squeezing operator 
followed by a displacement operator shifting the classi- 
cal position and momentum of the ion. The squeezing 
parameter is exactly the same as that of the centre-fixed 
squeezing operator in Sec. IIII CI for the same w 2 (£). A 
desired operation can be constructed by: first obtaining a 
t-o 2 (t) that produce the desired squeezing parameter; then 
finding a s(t) that can produce the desired displacement 
Pit) following the equation of motion Eq. (|A4|) . 



from the desired value g/\g\. The phase is possible to be 
rectified by first calculating the total phase shift gener- 
ated from our choice of w(t), then appropriate phase-shift 
operators are applied before and after the squeezing op- 
erator to compensate the extra phase shift. 

The condition of b(-T/2) = b(T/2) = 1 implies that 
uj 2 _ = luq before and after the operation. According to 
Eq. (|52p . this variation of u) 2 _ it) implies that the varying 
double well potential would bring two ions to proximity 
and then separated, as we have expected at the construc- 
tion. 



Appendix B: Construction of displacement operator 



The variation of the harmonic trap centre, s(t), that 
produces a specific displacement cto can be obtained sys- 
tematically by the inverse engineering method in Ref. 
[24| . We here present a simpler method that employs 
the linearity of displacements and paths. First of all, 
two arbitrary paths are speculated, Si(t) and S2(i); both 
paths satisfy the boundary conditions s(0) = s(T) = 0. 
According to Eq. (J25]), the paths will produce two dis- 
placements, «i = a[si] and a 2 = a^]. We require a\ 
and a 2 are not scaled by a real number, otherwise an- 
other path s^it) has to be speculated. If the requirement 
is satisfied, there must exist two parameters, 71 and 72, 
such that 



Appendix D: Construction of beam splitter 



a = 71 a 1 + 72« 2 



(Bl) 



Then the path s(t) = 7iSi(t) + 72^2^) will give the de- 
sired displacement ao- 



Appendix C: Construction of squeezing operator 

For a particular the auxiliary function bit) should 
behave as Eq. (fH|) at t > T after integrating Eq. (J7J with 
the boundary condition b(t < 0) = 1. To construct a 
squeezing operator with the desired squeezing operator 
g, we have to find an oj(t) that generates the parameter 
S = 2\g\ for bit > T). Such a condition can be satisfied 
by a wide range of w(i); a particular w(i) can be obtained 
inversely from a constructed b{t). An example is 

b s (t) = cosh 5 + sinhS sm(2ujt)h{t) + (I- h{t)) , (CI) 

where h(0) = and h(T) = 1. uj(t) should be continu- 
ous before and after the operation, so bsit), bsit), and 
bsit) have to be continuous at t = and t — T. For 
instance, h(t) = 10(t/T) 3 - 15(i/T) 4 + 6(i/T) 6 meets the 
requirement. The time variation of u>{t) can be obtained 
by inputting bs{t) into Eq. (J7J. 

So far we have neglected the phase of the squeeze 
operator. In the construction described above, the 
phase of the operation originates in the phase of (iT) 
and rjiT), as well as the phase shift in Eq. i.e. 
exp(i(0(T) + uj a T)). The u(t) is the unique factor de- 
termining the total phase, which is in general different 



Similar to the discussion in Sec. IIII Bl there are count- 
less forms of to- it) that can achieve a phase-shift operator 
on the breathing mode. It is generally difficult to guess 
the trapping potential magnitudes, A{£) and Bit), that 
can produce an appropriate uj-{t)\ they could be cho- 
sen only by trial-and-error. Alternatively, we outline the 
procedure of acquiring A{t) and Bit) inversely from an 
intellectually speculated bit). 

Consider the operation of double well runs from — T/2 
to T/2. bit < -T/2) and b(t > T/2) are the necessary 
condition for a phase-shift operator. An example would 
be 



b B {t) = 1 - fee 



-t 2 /o 



(Dl) 



where a <C T determines the time scale of the operation; 
the value k is chosen to generate the desired phase shift. 

The corresponding can be deduced inversely 

from Eq. ([?])• The time variation of the ion separation, 
r(t), is then obtained from i>j 2 _ (i) by Eq. (|52l) . A con- 
straint on Ait) and Bit) is obtained using the equation 
of motion Eq. (jT7)l and the time variation of r(t). To- 
gether with Eq. (jM)) . the unique form of A{t) and B(t) 
can then be found. 



Appendix E: Accuracy of beam splitter 

The quality of the beam splitter operation would be 
degraded by the anharmonic composition of the Hamil- 
tonian Hb, which are higher than second order quan- 
tum fluctuation that do not covered by the approximated 
Hamiltonian H 2 j i- e - 



Hb = H 2 - 



V2e 2 q 3 



2ire r 3 ir + V2g_) 



+ V2A(t)ri3q 2 + q- + q 3 _) 



Qf+e+qt) • 



(El) 



When comparing with H 2 , the magnitude of the anhar- 



monic terms is roughly offset by 



that scales as 



y/n, where n is the average number of phonon involved 
in the simulation. 
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process time is about only 20/o;o ~ 3.1 /is. In addition, 
the phase errors are less than 1% in all the runs. We 
conclude that for a bosonic simulation with single digit 
number of phonons in each mode, a quality phonon beam 
splitter can be implemented within a few /is. 

The accuracy of the beam splitter is worsened, as ex- 
pected, when more phonons are involved, because the 
factor */ (<7 2 ) /f increases. A higher operation speed also 
exacerbates the error. When considering Eq. (|E1[) in the 
interaction picture with respect to 7^2, the terms that are 
third order to */ (<7 2 ) /t are off-resonant. These nonlinear 
terms' contributions are suppressed by the RWA if the 
terms' variation is slow, then the effective Hamiltonian 
would be reduced to the forth order of yj (q 2 )/r. On the 
other hand, the off-resonant terms are significant for a 
high speed operation because the RWA is less effective. 



FIG. 4: Fidelity of the phonon state after a 50:50 beam split- 
ter operation with ujqo 2 = 2, 3, 4, 5, 7, 9. Total time for the 
double well to bring ions from and back to the pick-up position 
(step II to IV in Fig. gj) are about 11, 13, 15, 17, 20, 22 x l/u 
respectively. The minimum separation between ions is about 
1 lo in all the six runs. The dotted line shows a benchmark 
of 0.99 fidelity. 



The actual error produced crucially depends on the 
b(i) chosen. To demonstrate the feasibility of our beam 
splitter, we numerically simulate a 50:50 beam splitter 
using the be (t) in Eq. (|D1[) . Only the steps involving 
the double well are investigated (step II-IV in Fig. [3]); 
the pick-up process and the harmonic well transporta- 
tion are assumed to be error-free. The evolution of 
the motional states is tracked by numerically integrat- 
ing the Schrodinger equation with the Hamiltonian in 
Eq. (|E1[) . For simplicity, the interaction terms between 
the + and - modes are replaced by the expectation value, 
e.g. q\ —> (q 2 .), in the evolution of - mode; this treat- 
ment is accurate in our case because the back reaction 
scales as yj (q 2 )/r that is small. 

We consider the ions are 40 Ca + and the trapping fre- 
quency of the storage traps are wo = 27r MHz. The pick- 
up position is set as 50 Iq from the mid-point of the ions, 
where lo = y/ e 2 / 4ne muj 2 « 4.45 fxm is the ions' sepa- 
ration if they are placed in a single harmonic well. The 
operation speed is adjusted by tuning the characteristics 
time a, and the value of k is respectively chosen for each 
a to generate the desired phase shift. A quality beam 
splitter should produce a final state whose phase and fi- 
delity are close to those in the ideal case. 

We set the input states of both the centre-of-mass 
mode and the breathing mode to be pure Fock states, 
i.e. \tp(—T/2)) — |n+)|n_), such that the final states 
should be the same as the input states up to a phase. 
We find that the output fidelity of any |n + )|n_) input is 
higher than that of a | max{n + , max{n + , n_}} in- 

put, so only the runs with n + — n_ are shown in Fig. 0] 
for comparison. The double well process generates less 
than 1% error when uj 2 <j 2 > 7 for n max < 8; the total 



Appendix F: Ion separation without heating 

Usual processes of ion separation is to first reduce the 
trapping potential, then an additional potential is raised 
in the middle of the trap [37]. A fast process would 
heat up the ion significantly. In quantum computation, 
the side effect of heating can be compensated by subse- 
quent cooling that do not disturb the quantum informa- 
tion encoded in the electronic states. However, the time 
taken by ion cooling limits the operation speed of the ion 
trap quantum computer. Here we propose an alternative 
method to separate the ions, which can be rapidly con- 
ducted with minimal heating, by using an extra quartic 
potential. 

The configuration we are considering is the same as 
that in Sec. IIV1 where the Hamiltonian of the system is 
given in Eq. (|45[) . The ions are initially in the motional 
ground state at a common harmonic well with trapping 
frequency ujq, i.e. A(0) = and B(0) = towq/2. Our 
aim is to separate the ions to distant storage traps with 
frequency ujq with no motional excitation apart from the 
errors caused by the nonlinear terms. The initial centre- 
of-mass mode frequency is ujq and the breathing mode 
frequency is y/3uo 38]. Similar to the phonon beam 
splitter, we require 3A(t)r 2 /m + 2B(t) jm = lo 2 such that 
there is no excitation in the centre-of-mass mode. For the 
breathing mode, it can be shown that b(t < 0) = l/v^3 
is equivalent to the case of a single harmonic trap. Any 
b(t) with the boundary conditions 

b(t < 0) = -L , b(t > T) = 1, (Fl) 

can give the potential variation that retain the breath- 
ing mode in the ground state. b[t) and b(t) should also 
be continuous throughout the process unless r changes 
extremely fast. An example is 



3 -t ja 



(F2) 
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FIG. 5: Variations of potentials during heatingless ion separa- 
tion. Step I, a quartic potential is added to the common trap 
to form a double well. Step II, double well extends and sep- 
arates the ions to pick-up position. Step III, harmonic wells 
pick-up the ions and bring them to storage traps. 

where a is some characteristic time scale of the operation. 

After the ions' separation is large enough that the 
Coulomb interaction is negligible, the double well po- 
tential can be switched off and then the ions are picked 
up by two moving harmonic well. The velocities of the 
ions at the pick-up locations can be calculated by time- 
differentiating Eq. (|52p once with the oot(t) obtained 
from Eq. ([7]). The two displacing harmonic wells, which 
the trapping frequencies are ujq , should decelerate the ion 
motion to rest at the storage traps; the time variation of 
the potential can be deduced from inverse engineering 
method in Ref. [24j. The whole process is summarized 
in Fig. [5] Inverting the whole ion separation process can 
bring two individually trapped ions to a common trap. 



Similar to the phonon beam splitter, the error of this 
ion separation method is produced by the anharmonic- 
ity of the quartic potential and the Coulomb interaction, 
as shown in Eq. (jEll) . The magnitude of the generated 
errors highly depends on the b(t) chosen, but a faster 
ion separation process would generally create more error 
because the RWA becomes less effective. 

To demonstrate the feasibility of our scheme, the evolu- 
tion equation Eq. (|E1[) is numerically integrated using the 
&E(t) in Eq. (|F2|) with a — 2. For simplicity, we assume 
the operators of the other mode can be approximated by 
the expectation value, i.e. q± « (q±) in the equation of 
the =F mode. This assumption is accurate because the 
back reaction only produces higher order corrections to 
the to the already small errors. Our calculation simu- 
lates only the process of the double well extension (step 
I to II in Fig. [SJ heating effects caused by the pick-up 
and harmonic well transportation are negligible if the 
processes are well controlled. We consider the ions are 
Ca + and ojo = 2tt x 10 6 Hz, where the initial separation 
ro ~ 4.45 /im. After separating the ions from tq to about 
100 ro, only about 0.001 quanta is excited for both the 
+ mode and the - mode. The duration of this process is 
about 5 x 1/wq ~ 0.8/is. Because the transportation of 
the ions from the pick-up position to the storage traps 
(step III in Fig. [SJ can be arbitrarily fast [IH, the time 
duration of the whole ion separation process can be as 
short as 1 fis with minimal heating. 
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